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Abstract. We present general considerations regarding the 
derivation of the radial distances of coronal mass ejections 
(CMEs) from elongation angle measurements such as those 
provided by SECCHI and SMEI, focusing on measurements 
in the Heliospheric Imager 2 (HI-2) field of view (i.e. past 
0.3 AU). This study is based on a three-dimensional (3-D) 
magneto-hydrodynamics (MHD) simulation of two CMEs 
observed by SECCHI on Januai-y 24-27, 2007. Having a 3-D 
simulation with synthetic HI images, we are able to com- 
pare the two basic methods used to derive CME positions 
from elongation angles, the so-called "Point-P" and "Fixed- 
<j)" approximations. We confirm, following similar works, 
that both methods, while valid in the most inner heliosphere, 
yield increasingly large errors in HI-2 field of view for fast 
and wide CMEs. Using a simple model of a CME as an ex- 
panding self-similar sphere, we derive an analytical relation- 
ship between elongation angles and radial distances for wide 
CMEs. This relationship is simply the harmonic mean of the 
"Point-P" and "Fixed-0" approximations and it is aimed at 
complementing 3-D fitting of CMEs by cone models or flux 
rope shapes. It proves better at getting the kinematics of the 
simulated CME right when we compare the results of our 
line-of-sights to the MHD simulation. Based on this approx- 
imation, we re-analyze the J-maps (time-elongation maps) in 
January 26-27, 2007 and present the first observational evi- 
dence that the merging of CMEs is associated with a momen- 
tum exchange from the faster ejection to the slower one due 
to the propagation of the shock wave associated with the fast 
eruption through the slow eruption. 

Keywords. Interplanetary shocks (2139); Flares and mass 
ejections (7519); Instruments and techniques (7594) 



1 Motivation 

With the launches of the two Solar Terrestrial Relations 
Observatory (STEREO) spacecraft and the Coriolis space- 
craft in 2006 and 2003, respectively, coronal mass ejections 
(CMEs) can be, for the first-time, imaged continuously from 
the solar surface to 1 AU with coronagraphic and helio- 
spheric imagers. The CME on January 25, 2007 was the 
fastest eruption imaged by the STEREO/Sun-Earth Connec- 
tion Coronal and Heliospheric Investigation (SECCHI) suite 
to date (Howard et al. 20081). Although there was a 20-hr 



data gap in SECCHI coverage at the time of the ejection, 
these observations provide one of the best available tests for 
methods aimed at deriving CME dynamics from SECCHI 
observations for two main reasons. 

First, in contrast to slow ejections which arrive at Earth 
with speeds comparable to that of the ambient solar wind, 
a CME with initial speed greater than 1,300 km s^' should 
remain faster than the ambient solar wind in the entire HI-2 
field of view (FOV). Because most of the models of CME 
deceleration invoke a "drag" term proportional to the differ- 
ence between the ejection and the ambient solar wind speeds 



(Cargill 2004; Tappin 2006), the acceleration profile cannot 
be well constrained by the analysis of slow CMEs. 

A second reason is the presence of a preceding ejection 
from the same active region. This ejection was launched 16.5 
hours earlier and had a speed of about 600 km s"'^. Accord- 
ing to previous analyses ( [Lugaz et al. 



2009 Webb et al. 



|2009[ [Harrison et al.| |2069), the two eruptions interacted 
in the heliosphere somewhere between 20° and 30° elonga- 
tion from the Sun. It is expected that fast shock waves can 



propagate inside preceding ejections (Schmidt and Cargill 



2004[ Lugaz et al. 2005 1 and merge with the preceding shock 



waves. However, the variation of the shock speed inside the 
preceding magnetic cloud(s) is not known precisely. Nu- 
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merical simulations ( Lugaz et al. 2005 i have shown it can 
vary greatly due to the large variation in density, magnetic 
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field and alfvenic speed inside the magnetic cloud. There- 
fore, a constant or near constant speed cannot be assumed for 
the January 25, 2007 CME; in fact, most methods tested by 
[Webb et al.| ( |2009| l to explain the measurements, including 
cone models and numerical simulations, fared quite poorly 
past 25-30° (see their Figure 7) for at least one of the two 
observed fronts, although the cone model proved quite accu- 
rate in fitting the faster front. The authors noted that "con- 
version techniques from distance to elongation may require 
more work." It is the goal of this article to continue this pro- 
cess in an attempt to analyze HI observations better. 



2 SECCHI observations of the January 24-25 CMEs 
and numerical simulation 

The two successive CMEs of January 24-25, 2007 were ini- 



tially reported by Harrison et al. (2008 1. At the time, the two 
STEREO spacecraft were still in close proximity with Earth 
(within 0.5°) and STEREO-A was rolled by about 22° from 
solar north. Beyond COR-2 FOV, only STEREO-A/SECCHI 
observed these eruptions originating from an active region 
behind the eastern limb. The two eruptions were first imaged 
by COR-1 at 14:03UT on Januaiy 24, 2007 and 06:43UT 
on January 25, 2007. Based on their appearance in corona- 



graphic images, we determined in Lugaz et al. ( 2009 1 that 
they were associated with active region 10940 which was 
about 20° behind the eastern limb at the time of the first erup- 
tion. Due to positions of the Solar and Heliospheric Obser- 
vatory (SOHO) and STEREO spacecraft in January 2007, no 
triangulation of the source region of the eruptions is possible, 
as was done for later CMEs by Howard and Tappin ( 2008| l for 
example. 

Based on the time-height profiles of the CMEs in the 
SOHO/LASCO FOV, and using the same position angle (PA 
90) for both CMEs, the speed in the corona of the first CME 
was 600 km s^^ and it was 1350 km s^^ for the second one. 
The data gap in SECCHI coverage started after 04:53UT 
and 09:53UT on January 25, 2007 for STEREO-A and B, 
respectively and lasted until the start of January 26, 2007. 
Assuming no deceleration, the two ejections should have in- 
teracted during this time. After the SECCHI data gap, two or 
three bright fronts associated with the eruptions were tracked 
in in-2 (Harris on et~aLl [2008] [Lugaz et alf |200'8| |2009t 
Webb et al.^ g009 1, the first front up to elongation angles of 



about 55° with HI-2 and up to much larger elongation angles 
(- 90°) with SMEI ( Webb et al.|[2009] l. SMEI observations 
could not help during the SECCHI downtime, because the 
CMEs were inside the SMEI exclusion zone circle of 20° 
around the Sun. 

We performed a numerical simulation of these ejections 
with the Space Weather Modeling Framework (SWMF) 
(|T6th et al. 2005) using the solar wind model of jCohenj 
|et al.| ( (2007^ . The simulation set-up and detailed results have 
been pubUshed inlLugaz et al. I (1200911. Simulations with the 



ENLIL model of |Odstrcil et"aL] ( |2005l ) and the HAFv.2 model 
of Hakamada and Akasofu ( 1982[ i and Fry et al. ( 200 1[) hav e 
also been performed and published in Webb et al. ( 2009 1. 
Based on the numerical analyses, the fronts observed by HI- 
2 and SMEI have been associated with the two CMEs, val- 
idating the numerical models on one hand and helping the 
analysis of the complex observations on the other hand. The 
goal of the current study is to test the existing methods to 
derive CME radial distances from elongation angles with the 
help of a 3-D simulation. 



3 Determining CME positions from elongation angles: 
Testing the existing methods 

So far, CME positions have been determined from STEREO 
observations via 3-D forward fitting of a cone-model or a 



flux-rope- shaped density enhancement ( Boursier et al. 2009 



[Thernisien et aH [2009) , via 3-D reconstruction in COR-2 
FOV (i.e. within 20 Rq) ( [Mierla et al.] [2008] [de Koning 
et al.1'2009), by mass conservation principles ("Colaninno and 



Vourlidas 2009 ), or by applying one of two simple approxi- 
mations giving an analytical relation between elongation an- 



gles and CME positions ( Wood et al. 2009 Rouillard et al. 



[2009[ [Davis et al.|[2009[ ). These analytical relations provide 
a quick and easy way to estimate CME dynamics in the he- 
liosphere. 3-D reconstruction and forward modeling are ex- 
pected to be more accurate than these simple relations, espe- 
cially in the COR FOV where they have been mostly used so 
far, but they also have some limitations. For example, the 3-D 
reconstruction methods require multiple viewpoints, which 
might become less and less frequent as the STEREO space- 
craft separate; when there are multiple observations, they as- 
sume that both SECCHI instruments observe the same struc- 
ture, which is not true in the HI FOV. Additionally, forward 
modeling attempts to fit geometrical and kinematic informa- 
tion at the same time. To simplify the fit, a kinematic model 
(often constant acceleration or constant speed) is usually as- 
sumed. As noted above, these assumptions cannot be used 
for complex events, such as those involving CME-CME in- 
teractions. 

3.1 The "Point-P" and "Fixed-0" approximations 

The intensity of the Thomson scattering depends on the an- 
gle between the scattering electron, the Sun and the observer 
(Mi nnaertj|1930| l. The loci of the ensemble of points where 
the intensity of Thomson scattered light is maximum is re- 
ferred to as the "Thomson surface" ( [Vourlidas and Howard| 
2006). In 3-D space, this surface lies on the surface of a 
sphere with the Sun-observer line as the diameter, and so we 
refer to this as the Thomson sphere from now on. A simple 
plane-of-the-sky approximation cannot be used with accu- 



racy in the HI FOV (e.g., see VourUdas and Howard 20061. 
Therefore, to know which part of a CME is imaged, one 
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Fig. 1. Left: Geometry of the observations and the methods described in the article. The figure corresponds to the January 24-25 CMEs in 
the late phase of their merging. This illustrates the different CME positions obtained from one measurement of the elongation angle at e. 
The black and yellow circles illustrate the model of CMEs used to derive the relation described in the article and the Point-P approximation, 
respectively; the white circle is the Thomson sphere; the green dot and black disk are STEREO-A and the Sun, respectively (not to scale). 
The angle (j> is set at 90° to determine the CME distances but it is shown here as determined from the position of the active region at the 
start of the eruptions. Right: The model (expanding propagating sphere) proposed to derive CME position from SECCHI measurements is 
illustrated for a different simulation (August 24, 2002) with different models of the solar wind and CME initiation. The yellow sphere is 
centered at the Sun, the white translucent sphere is the model of the CME front and the actual simulated CME is shown as an isosurface of 
scaled density 20 cm~'^ AU^'^ color-coded with the speed. 



needs to consider the complex interaction of the CME 3-D 



density structure with the Thomson sphere (e.g., see Lugaz 
|et al.[[2008| . An additional problem is that the speed and ac- 
celeration should be calculated for the same plasma element 
(i.e. usually for a single radial trajectory). Even if the CME 
positions can be determined accurately from HI observations, 
further assumptions regarding the CME geometry must be 
made to derive kinematic information, since what is observed 
over time is not necessarily the same part of the CME, as 
shown in [Lugaz et al.] ( p009] l and [Webb et al.| ("2009 ) . There 
are two main simple approximations which have been used 
to replace the plane-of-sky approximation for heliospheric 
measurements: they are referred as "Point-P" and "Fixed-t/)" 
(iKahler and Webbl [2007] [Howard et al.| [20071 |Wood eFaLj 
2009| l; the geometry of the observations and the reconstruc- 
tion is illustrated in the left panel of Figure 1 for a plot of the 
simulations of the January 24-25, 2007 CMEs. 

The "Point-P" (PP) approximation is the simplest possi- 
ble way to relate elongation angles to CME radial distances 
while taking into account the Thomson sphere geometry. As- 
suming a spherical front centered at the Sun, the elongation 
angle e and the position of the CME i?pp are related by: 



where c?stereo i'^ 0.97 AU) is the heliocentric distance of 
STEREO-A for this event. The CME front obtained from 
this approximation is shown with the yellow circle in the left 
panel of Figure 1. Obviously, this approximation is poor 



for narrow CMEs such as the one studied by Wood et al. 



( 2009 1 and for dense streams and corotating interacting re- 



i?PP 



^STEREO sm e, 



gions (CIRs) which are structures of narrow azimuthal ex- 
tent at 1 AU (the typical width is less than 20° as inferred 
from |Jian et al.[ ( [2006| l for example). Even for wide CMEs, 
the CME fronts are not spherically symmetric, in part due 
to their interaction with the structured coronal magnetic field 
and solar wind. This has been shown by multiple-spacecraft 
observations (e.g. 'Mostl et al. 2009) and from simulations 
([Riley et al.[ [2003, Manchester et al., 2004, Odstrcil et aT] 
2005| l. Last but not least, the reconstructed CME position is 
independent of the propagation angle! Webb et aL] ( [2009[ l re- 
marked that the PP approximation is not adequate far from 
the Sun, e.g., in HI-2 and SMEI/camera 2 FOVs. 

The "Fixed-(/)" (F0) approximation, in turn, takes the op- 
posite philosophy and considers that a single particle, prop- 
agating on a fixed-radial trajectory, is responsible for the 
Thomson scattered light. The elongation angle measurement 
must simply be "de-projected" from the Thomson sphere 
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onto this radial trajectory, resulting in the relation 



200 - 



R 



F4> 



ds 



sine 



TEREO 



sin(e + ( 



where ip is the angle between the Sun-observer line and the 
trajectory of the particle. The position obtained from this ap- 
proximation is noted as Rp^f, in the left panel of Figure 1 . Ob- 



viously, this approximation is well adapted for CIRs (Rouil 



lard et al. 2008) and small "blobs" ( [Sheeley et al.[ |2008 



2009(,Rouillard et al.| 12009) . However, since it assumes that 
what is tracked is a single point, the method is expected to 
work poorly for wide CMEs. The equation can be fitted for 
(p (assuming no or constant acceleration), giving the origin 
of the transient ( Rouillard et al.[ 2009| l and/or the speed. The 
main limitation of this method is that it completely ignores 
the CME geometry. It also does not take into account the 
angle dependency of the Thomson scattering. 

3.2 Comparison with 3-D simulated data 

We test the two methods with our synthetic line-of-sight pro- 
cedure and compare the resulting positions to the 3-D sim- 
ulation for the second CME (January 25 CME) at PA 90. 



This work is the continuation of section 4.3 from Lugaz et al. 
( |2009) l and its associated Figure 6. We derive the elongation 
angles and radial distances of the CME front as follows: for 
the line-of-sight images, we use elongation angles measured 
at the point of maximum brightness at PA 90. For the nu- 
merical simulation, we use the position of maximum density 
along different radial trajectories (at longitudes 90°, 80° and 
70° east of the Sun-Earth line) and on the Thomson sphere, 
all of these in the ecliptic plane (PA 90). Results are shown 
in the top panel of Figure 2. 

Below approximately 100 R©, the two methods give sim- 
ilar results differing by less than 10%. The Fcf) method 
gives slightly better results when compared to the nose of 
the CME; the PP approximation works best if one assumes 
it tracks the intersection of the CME front with the Thom- 
son sphere (see middle panel of Figure 2). Above 100 Rq, 
the two methods give increasingly different results. Com- 
pared to the simulation results along all three radial trajec- 
tories presented here, the PP approximation results in a too 
large deceleration of the CME, whereas the F0 results in an 
apparent acceleration. This acceleration appears unphysical, 
since CMEs faster than the ambient solar wind are expected 
to monotonously decelerate due to a "drag" force. Simi- 
lar results have been reported, most recently by [Wood et al.| 
2009| l. The middle panel of Figure 2 shows the errors be- 



tween the position of the CME front at the limb and the po- 
sition from each of the two methods, as well as the error be- 
tween the PP position and the intersection of the CME front 
and the Thomson sphere. Although the errors are fairly low, 
they can result in large errors in the velocity and acceleration 
of the CME (see bottom panel of Figure 2). These methods 
can provide an average speed of the CME front within the 




Time (hr) 



Fig. 2. Position (top), error (center) and speed (bottom) of the sec- 
ond CME front at PA 90 from the simulation and as derived from 
the synthetic SECCHI images with the different methods. 



first 100 Rq, but they cannot be relied upon to study com- 
plex physical mechanisms such as CME -CME interaction. 
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4 Improved method to determine CME position 

Based on the relatively poor results for the PP and F(j> meth- 
ods, we propose another analytical method based on simple 
geometric considerations and a simple model of CMEs. We 
construct this model on a few principles: first, it should take 
into account the geometry associated with the Thomson scat- 
tering as well as the CME propagation and second, it should 
have the lowest number of free parameters possible. To con- 
struct such a model, we start from the knowledge that CMEs 
are known to evolve self-similarly in the heliosphere (e.g., 
see Krall et al. 2006 1. The simplest approximation is to as- 



sume that the CME peak density maps out as a sphere con- 
nected to the center of the Sun; the center of the sphere prop- 
agates in a fixed, radial trajectory (see right panel of Fig- 
ure 1). In contrast to the PP approximation, the sphere is 
not centered on the Sun. Consequently, this method takes 
into account the direction of propagation of the CME. This 
approximation is also the one used in |Webb et aL] ( |2009| l to 
produce their Figure lb. 

If we assume no deflection of the CME in the corona or 
in the heliosphere, the angles defining the trajectory of the 
center of the sphere can be derived from the flare informa- 
tion (with an understanding of the limitations in the connec- 
tion between flares and CMEs) or from forward modeling 
of the COR observations or mass analysis. For the January 
25, 2007 CMEs, we will consider that the center propagates 
from the eastern limb at PA 90. There are many ways this 
sphere "interacts" with the Thomson sphere to produce the 
Thomson-scattered signal. We consider two hypotheses: the 
geometry associated with the Thomson scattering is dom- 
inant and the emission originates from the intersection of 
the sphere with the Thomson sphere or it is negligible and 
the emission originates from the line-of-sight tangent to the 
sphere (see left panel of Figure 1 for the geometry and the 
notation used). The first hypothesis gives di — Rp^ for the 
diameter of the circle representing the CME front at the PA 
where the measurement is made. This PA can differ from 
the latitude A along which the center of the CME propagates. 
Correcting for this, the nose of the CME is at a distance of 
i?F(^/cos(PA — A). This gives a new interpretation for the 
"Fixed-(/)" approximation, namely that it gives the diameter 
of the circle representing the CME at each PA, assuming the 
emission originates from the intersection of this circle with 
the Thomson sphere. 

The distance of the point tangent to the CME along the 
given PA (see the left panel of Figure 1 for the notation) is: 



dsin e 



cos a 



dsin e 



f)) 



The diameter of the circle representing the CME at this PA is 
simply given by: 



dsin e 



which is the harmonic mean of the PP and F0 approxima- 
tions. To obtain the diameter of the sphere, this must also 
be corrected for the difference between the measured PA and 
the direction of propagation of the CME: 



1 _ cos(PA - A) / 1 
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1 



R 



F0 



R 



PP 



COS a 1 + sin(e - 



This correction is required because all parts of a CME cannot 
be assumed to move radially outward with the same speed. 
Thus, this hypothesis is most likely true for the nose of the 
CME, which is where the speed must be calculated. We plot 
the position, error and speed derived from this approximation 
(referred as the harmonic mean (HM) approximation) in the 
three panels of Figure 2. As can be seen, this simple model 
gives better results than the PP and Fcf) approximations, espe- 
cially for the speed of the CME at large elongation angles. 



5 Revisiting tiie January 26-27, 2007 observations: 
CME-CME merging 

5.1 Data analysis 

With this method, we re-analyze the data from the two fronts 
observed by HI-2 on January 26-27, 2007. We analyze the 
data at PA 69, where the SECCHI's coverage is best for this 
event. There were only limited observations of the two ejec- 
tions prior to the data gap. For the first ejection, all three 
methods agree and give an average speed between 550 and 
600 km s^i at 40 R©, which is consistent with LASCO ob- 
servations and also with the speed of 604 km s^^ reported by 
[Harrison et al.| ( |2008| l for the front at PA 90 (i.e. the nose of 
the ejection). For the second eruption, we use LASCO data, 
which give a speed of approximatively 1,200-1,300 km s^^ 
at 20 R0. Next, we analyze the two fronts after the data gap 
in HI-2 FOV. The top panel of Figure 3 shows the derived po- 
sition for the two fronts according to the three methods. First, 
it is worth noting that the Fcf) and the HM approximations dif- 
fer by less than 10% up to approximatively 180 Rq (40°), but 
the HM approximation does not result in a large apparent ac- 
celeration at very large elongation angles. Next, we derive 
the speed of the two fronts according to these methods. We 
plot a running average over approximatively 5 hours to re- 
duce the magnitude of the error in the speed. HI-2 resolution 
is 4 arcmin; assuming the elongation angles are measured 
with a precision of 5 pixels, the error in position is of the 
order of 1.5% and the resulting speed has an error of about 
15%. 

We beUeve the analysis of the numerical results from sec- 
tion 3.2 shows that the PP method cannot be used to study 
the speed of limb CMEs past 100 Rq, which is the approxi- 
mate position of the two fronts after the data gap. According 
to the Fcf) and HM methods, the second front, which fades 
out at about 33° elongation, has an average speed of 680 
and 605 km s"^, respectively, with a general decelerating 
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Table 1. Summary of the speeds measured by SECCHI for the two 
fronts. 
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Fig. 3. Position (top) and speed (middle and bottom) of the two 
fronts at PA 69 according to the three methods. The errors are typi- 
cally 1.5% for the position and 15% for the speed. The averages are 
shown with dotted lines and the second front with dashed lines. 



trend with an initial speed around 750-850 km s~^ around 
100 Rq. The two methods are overall consistent with each 
other, and we believe this shows that the transient associated 
with the second front had a speed of 750-850 km s^^ around 
100 R0 and decelerated to 500-600 km s^^ before disappear- 
ing around 140 R0. 

For the first front, which is tracked until 53°, the F(j) re- 



Front 


Speed Before Collision 


Speed After Collision 


1 


600 km s"^ 


850-900 km s"^ 


2 


1200-1300 km s"^ 


800 km s"^ @ 80 R© 






550kms"^@ 140 R© 



suits in a strong acceleration after 40° elongation and the PP 
method in an almost constant low speed. Fast CMEs are not 
expected to experience large acceleration in the heliosphere 

|200§. The HM 



(e.g. Gopalswamy et al. 2001 Tappin 



method results in a speed more consistent with this fact than 
the F(j) method, although it shows a limited, unphysical ac- 
celeration at large elongation angles. The average speed ob- 
tained from the three methods is 490, 1340 and 845 km s^^, 
respectively; the average speed of the F(p and HM methods 
for observations between 28° and 40° is 880 km s^^ and 
705 km s~^. The analysis is more complicated than for the 
second front, but, we believe that the observations are con- 
sistent with a transient whose average speed is about 850- 
900 km s^^ (the average value of the HM method, and the 
average value of the F(/) within 40°). 

5.2 Consequence for the process of CME-CME interac- 
tion 

The derived speeds of the fronts are summarized in Table 1 . 
We believe there are 4 scenarios consistent with the result 
that the first front is faster than the second front after the data 
gap; we analyze these scenarios with respect to the measured 
speeds of the two fronts. A schematic view of the 4 possibili- 
ties is shown in Figure 4. In the first scenario, the January 25 
CME could have "passed" the January 24 CME without ma- 
jor interaction. This scenario is possible if the two eruptions 
have a large angular separation, and if they do not propagate 
along the same direction. Then, part of the fast front could, in 
the projected images, "pass" the slow front when in fact there 
is no interaction. This scenario is described in greater details 
in |Webb et al.| ( |2009| ). While it is plausible that only a small 
part of the two CMEs interacted and that the major part of the 
January 25 CME simply passed next to the January 24 CME 
without interaction, we believe this is very unlikely. First, it 
is hard to understand how the speed of the January 24 erup- 
tion could be faster after the data gap than before; also, the 
January 25 eruption shows a strong deceleration during the 
data gap, which tends to suggest some form of interaction. 
Second, the measured width of the eruptions -greater than 
100° in LASCO FOV as reported in |Webb et al.] ( |2009l l- also 
makes a missed encounter implausible. Last, this is not sup- 
ported by any MHD models, which tend to show that CMEs 
act as magnetic barriers. This scenario could however ex- 
plain what happened if the two CMEs were associated with 
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Fig. 4. The four scenarios for CME-CME collisions that might explain the fact that the first front after the collision is faster than the second 
front. In the sketches, the ellipses, the solid arcs, and the dashed arcs correspond to the ejecta, dense sheaths and to the shock waves, 
respectively. 



different active regions and, consequently, had a large angu- 
lar separation. This separation could be as large as 35° if 
the first CME was associated with the eastern most active 
region and the second CME with the western most active re- 
gion present in January 24-25, 2007. Our arguments to asso- 
ciate both ejections with the same active region can be found 
inP 



Lugaz et al. (2009 1. 



In the second and third scenario, the two CMEs collide, 
the collision is associated with momentum transfer between 
the ejections (as Farmgia and Berdichevsky ] 12004 ) consid- 
ered). The observations appear to be consistent with both 
eruptions having the same speed after the collision, i.e. a 
perfectly inelastic collision. However, it is hard to under- 
stand the evolution of the speed of the two CMEs after the 
collision according to this scenario. If the January 25 CME 
pushes the January 24 CME, both fronts should have a sim- 
ilar speed at all times after the collision. This scenario ap- 
pears more plausible if one believes the speeds derived using 
the PP method. However, using the PP speeds and positions, 
the average transit speed of the two fronts during the data 
gap should be 500 and 650 km s^^ respectively. This sce- 
nario would therefore be consistent with a large deceleration 



of the January 25 (fast) CME and almost no acceleration of 
the January 24 (slow) CME, which, in turn, can only hap- 
pen if the January 24 CME is much more massive than the 



January 25 CME. Webb et aL]( |2009| ) reported the mass of the 
January 24 and 25 CMEs being 4.3 x lO^^ g and 1.6 x 10^^ g, 
respectively, making this scenario very unlikely. 

In the third scenario, the collision is elastic and there is a 
momentum transfer from the second to the first ejection on 
a time-scale of 12-20 hours. The momentum transfer has an 
unknown cause and continues until the second eruption be- 
comes slower than the first one. This scenario is not funda- 
mentally different from the last one, which does not require 
unknown processes and can explain the disappearance of the 
second front. 

In the fourth scenario, the unknown process is, in fact, 
the compression and momentum transer associated with the 
shock wave from the January 25 CME. Before the CMEs col- 
lide, the shock wave driven by the January 25 CME propa- 
gates through the January 24 CME (ejecta and sheath), com- 
pressing and accelerating it, before merging with its associ- 
ated shock wave. After the data gap, the first front corre- 
sponds to the sheath associated with the merged shocks. Due 
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to its interaction with the January 24 CME and sheath, the 
shock wave initially associated with the January 25 CME has 
decelerated rapidly to a speed ^ 850 km s^^. There are two 
possibilities to explain the second front: it could be the rem- 
nant of the sheath associated with the January 25 shock wave 
which is "trapped" between the two CMEs and "forced" to 
propagate with a speed comparable to that of the January 24 
CME. However, the distance between the two fronts is of the 
order of 20 Rq at PA 69. If the January 24 CME is between 
the two fronts, this would mean that the magnetic cloud has 
been compressed to less than 20 Rq, which does not seem 
reasonable. Moreover, the two fronts appear to merge along 
PA 90, which is inconsistent with this explanation. The 
other possibility is that it is associated with a transient phe- 
nomenon during the shock-CME or shock-sheath interaction 
or three-dimensional effects. For example, the core (or any 
part of the cloud) of the January 24 CME could be have been 
compressed and accelerated by the shock wave and it relaxes 
to slower speeds. Most likely, part of the sheath associated 
with the January 24 CME gets compressed to very high den- 
sity and relaxes to the average value of the new sheath (sim- 
ilar to what has been discussed in|Lugaz et al. [2005 1). How- 
ever, each of these sub-scenarios involve the propagation of 
the January 25 shock through the January 24 CME. We note 
that this scenario does not require the presence of a shock 
wave driven by the January 24 CME, but simply a sheath of 
dense material (piled-up mass and/or compressed material) 
ahead of the CME. The only difference due to the possible 
absence of the first shock wave is that there is no instance of 
shock-shock merging. Therefore, the shock wave ahead of 
the merged CMEs after the interaction is simply the shock 
wave originally driven by the January 25 CME now propa- 
gating into an unperturbed solar wind. 



6 Discussions and Conclusions 

In the first part of this study, we have tested the two most 
common methods used to derive CME radial distances from 
elongation angle measurements, the Point-P and Fixed-cj) 



larger elongation angles. The Point-P and Fixed-(/) approx- 
imations are expected to yield a lower and upper-bound to 



methods. Confirming previous work by Kahler and Webb 
( |2007] l , |Wood et al.| ( |2U09] l and |Webb et al., (20091 ) we find 
that, above 35°, both methods yield poor results, especially 
for CME speed and acceleration. We propose an alterna- 
tive analytical method to derive CME radial distances. We 
consider a very simple model, namely that the density peak 
maps out as a sphere whose center propagates radially out- 
ward from the flare location, and that the elongation angle 
corresponds to the angle of the line-of-sight tangent to this 
sphere. We find that the diameter of this sphere is given by 
the Iiarmonic average of the Point-P and Fixed- </) approx- 
imations further corrected by l/cos(PAapp) where PAapp is 
the position angle with respect to the nose of the CME. For a 
limb ejection, this method gives results similar to the Fixed- 
(j) approximation up to about 40° and more realistic results at 



the actual distance of a CME (e.g., see Webb et al. 2009 1. 
Any alternative methods to determine radial distances from 
elongation angles shall fall in-between, as is the case here. 
However, we find a particular physical interpretation for the 
harmonic mean of these two methods. We have also found a 
new interpretation of the position derived from the Fixed-(f> 
approximation, namely that it is the diameter of the sphere 
representing the CME if the emission is assumed to originate 
directly from the Thomson sphere. This might explain why 
this approximation works fairly well even for wide CMEs. 

We must be aware of the limitations of this method. First, 
this is most appropriate for wide CMEs such as the ones ob- 
served in January 24-25, 2007, because the assumed geome- 
try is consistent with a CME whose angular extent is 90°. 
This approximation, while arbitrary, is required to reduce 
the number of free parameters of the model to one. It also 
appears to be a better approximation for wide CMEs than 
the Point-P approximation which is consistent with a CME 
whose angular extent is 360°. Secondly, this model assumes 
that the CME propagates on a fixed radial trajectory, ignoring 
heliospheric deflection. This is the same assumption made 
to derive the Fixed-i/) approximation and is also required to 
reduce the number of free parameters. In future works, we 
shall investigate how stereoscopic observations of CMEs by 
the two STEREO spacecraft can help relax these two condi- 
tions. Last, the model of CMEs use to derive this approxi- 
mation assumes that the CME front (part piled-up mass, part 
shocked material) maps out as a sphere. As noted in the in- 
troduction, CME fronts are known to be distorted and usu- 
ally flattened from their interaction with the structured solar 
wind. In Figure 1, we have shown two simulated instances 
where this approximation is more or less appropriate; it is 
worth noting that the two simulations use different models of 
CME initiation and solar wind. Assuming a more complex 
shape (for example a "pancake") would require a fitting of 
the model and could not yield a direct analytical relationship 
such as the one derived here. 

We have re-analyzed the HI-2 measurements in January 
26-27, 2007 associated with two interacting CMEs with 
the three methods. We found that the first bright front 
after the interaction corresponds to a transient propagat- 
ing with a speed of about 850 km s^^, while the second 
front corresponds to a transient whose speed decreases from 
850 km s^^ to 550 km s~^ in about 12-18 hours before ulti- 
mately disappearing. Among the 4 scenarios which could ex- 
plain that the acceleration of the first front relative to the sec- 
ond one, we found that the only scenario consistent with the 
observations require that the part of the shock wave driven 
by the (faster) January 25 CME propagates first through the 
(slower) January 24 CME. The propagation of this fast shock 
wave inside the CME and dense material of the sheath results 
in its large deceleration. The most likely explanation for the 
origin of the second front is that it is part of the dense sheath 
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associated with the January 24 ejection, which is compressed 
and accelerated by the shock and decelerates to the speed of 
the January 24 magnetic cloud. 

Our analysis has been limited to one case of a fast, wide 
limb CME. It is for this particular geometry that the Fixed-(^ 
approximation is expected to give the largest error at large 
elongation angles. However, we believe our new average 
method should provide an improvement over the existing 
methods, notably over the Point-P approximation for wide 
eruptions, and that it should be used compUmentary to three- 
dimensional fitting methods and numerical simulations. We 
plan to test and validate this relation for other heliospheric 
observations of wide and fast CMEs, starting in the near fu- 
ture with the April 26, 2008 eruption. Our analysis of the 
January 24-27, 2007 observations is the first hehospheric ob- 
servational evidence of a shock wave propagating inside a 
CME. More observations without data gaps are required be- 
fore we have a more definite understanding of CME-CME 
interaction. 
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